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DIGITAL EQUALISING METHOD 
AND RADIO COMMUNICATION RECEIVER IMPLEMENTING SAID METHOD 

The invention concerns the digital equalization of 
signals. It has important applications in the field of 
radio conmunications . 

The method is applied to the processing of a signal 
received from a transmitter over a channel between 
transmitter and receiver, whose response is known or has 
been estimated beforehand. One of the main issues is then 
that of the compromise between the performance of the 
equalizer and its complexity. 

A complete maximum likelihood estimation of the 
transmitted signal is possible, for example by employing 
the Viterbi algorithm (see G.D. Forney Jr. : The Viterbi 
Algorithm, Proc . of the IEEE, Vol. 61, No. 3, March 1973, 
pages 268-278) . Nevertheless, as soon as the impulse 
response of the channels becomes long or the number of 
possible discrete values becomes large, the exponentially 
increasing complexity of these methods renders them 
impracticable . 

We consider the case of a radio communications 
channel used for the transmission of a signal composed of 
successive sequences or frames of n symbols, dj^ (l<k<n) . 

The symbols dj^ have discrete values : binary (±1) in the 
case of a BPSK (binary phase shift keying) modulation; 
quaternary (±l±j ) in the case of a QPSK (quaternary phase 
shift keying) modulation. 

After baseband conversion, digitization and matched 
filtering, a vector Y of the received signal, 
corresponding to the symbols transmitted over the duration 
of a frame, is defined by the expression : 
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where W+1 is the length in numbers of bits of the 
estimated impulse response of the channel, 

r = (rQ, r^, ..•/ r^) is the estimated impulse response of 
the channel, the r^ being complex numbers such that 37^=0 
for q<0 or q>W, y^^ is the k-th complex sconple received 
with l<k<L=n+W, and Yj^ is a vector of dimension L composed 
of samples of additive noise y^^k* "^^^ estimated impulse 
response r takes into account the propagation channel, the 
signal shaping by the transmitter and the receiver 
filtering. 

The matrix A of dimension Lxn has a Toeplitz- type 
structure, meaning that if j denotes the term in the 
i-th row and in the j-th coliimn of the matrix A, then 
^i+l,j+l = l<i<L-l and l<i<n-l. The terms of the 

matrix A are given by : ajL^j=0 for l<j<n (A therefore has 
only zeros above its principal diagonal) ; fo^: 
W+l<i<L (band-matrix structure); and for l<i<W+l. 

The matrix equation (1) expresses the fact that 
the signal received Y is an observation, with an additive 
noise, of the convolution product between the channel 
impulse response and the transmitted symbols. This 
convolution product can also be expressed by its 
Z- trans form : 



Y(Z) = R(Z) ,D(Z) + Yjj(Z) (2) 

where D(Z) , Y(Z), R(Z) and Yj^{Z) are the Z-transforms of 

the transmitted symbols, of the received signal, of the 
impulse response and of the noise, respectively : 

n 

D(Z) = 2 ^k- (3) 
L 

Y(Z) = 2 yj,. Z'^ (4) 
k=l 

w 

R(Z) = X ^q- 2""^ (5) 
q=0 

A conventional solution to solve a system such as 
(1) is the so-called "zero forcing" method, by which we 
determine the vector D^p of n continuous components which 

minimizes the quadratic error 8 = ||ad-y||'^ , Subsequently, a 
discretization of the components of vector Dgp relating to 
each channel is performed, often through a channel 
decoder. The least- squares solution is given by : 

DgF^ (A^A) *^A^Y, where A^ denotes the conjugate transpose of 
matrix A. We then are faced with the problem of inverting 
the Hermitian positive matrix A^A. The inversion can be 
effected by various classical algorithms, either directly 
(method of Gauss, Cholesky, etc.) or by an iterative 
technique (Gauss -Seidel algorithm, gradient algorithm, 
etc . ) . 

The estimation error D-D^p is equal to (A^A) "^A^Yj^, 
hence the solution includes noise with a variance: 

= E(||D-D2Kr) = Trace[(A^A)"^j (6) 

where Nq is the noise power spectral density. It can be 



seen that noise enhancement occurs when the matrix A^A is 
badly conditioned, i.e. when it has one or more 
eigenvalues close to zero. 

This noise enhancement is the main drawback of the 
conventional solution methods. In practice, the cases 
where the matrix A^A is badly conditioned are frequent, 
especially with multiple propagation paths. 

A relatively simple means to partly remedy this 
drawback is known, by accepting residual interference in 
the solution, i.e. by not adopting the optimism least- 
squares solution, but the solution : Dj^gj,= (A%+ Nq ) "^A^, 
where Nq denotes an estimation of the noise spectral 
density, that the receiver must then calculate. This is 
known as the MMSE (minimum mean square error) method and 
it allows the estimation variance to be reduced relative 
to the zero- forcing method, but introduces a bias. 

The zero forcing methods and the like amount to 
performing an inverse filtering on the signal received by 
a filter which models the transfer function 1/R(Z) 
calculated using a certain approximation (quadratic in the 
case of zero forcing) . Where one or more roots of the 
polynomial R(Z) (equation (5)) are situated on the unit 
circle, the theoretical inverse filter presents 
singularities rendering it impossible to estimate by a 
satisfactory approximation. In the case of the quadratic 
approximation, this corresponds to a divergence in the 
error variance when the matrix A^A has an eigenvalue 

equal to zero (relation (6) ) . 

This problem does not arise in methods such as the 
Viterbi algorithm which intrinsically take the discrete 
nature of the symbols into account, but which require a 
much higher calculating power for large- sized systems. 
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U.S. Patent 4,701,936 discloses a channel equalizer 
using an all -pass filter determined with reference to the 
Z- transform of the estimated channel impulse response. 

An object of the present invention is to propose an 
ec[ualization method achieving a good compromise between 
the reliability of the estimations and the complexity of 
the equalizer . 

Another object is to produce an equalizer requiring 
a reasonable calculating power but capable of processing, 
^ with performance comparable to that of Viterbi equalizers, 
signals whose symbols have a relatively high nximber of 
states and/or signals carried on a channel with a 
relatively long impulse response. 

Accordingly, the invention proposes a digital 
eq[ualization method for estimating discrete information 
symbols from digital samples of a signal received over a 
transmission channel represented by a finite impulse 
response of W+1 coefficients, W being an integer greater 
than 1. This method comprises the steps of: 

- determining the W roots in the complex plane of the 
Z- transform of the impulse response; 

- distributing the W roots into a first set of W-p 
roots and a second set of p roots, p being an 
integer greater than 0 and smaller than W, the roots 
of the second set being closer to the unit circle 
than those of the first set according to a 
determined distance criterion in the complex plane; 

- obtaining an intermediate signal by applying a first 
equalization method to the received signal based on 
a finite impulse response whose Z-transform, formed 
by a polynomial of degree W-p in Z"^, has roots which 
are the W-p roots of the first set; and 

- obtaining estimations of the discrete information 
symbols by applying a second equalization method to 
the intermediate signal based on a finite 
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impulse response whose Z- transform, formed by a 
polynomial of degree p in Z'^, has roots which are 
the p roots of the second set. 

The "first equalisation method" will usually be 
chosen so as to treat the unknown symbols as continuous 
variables. This leads to a process similar to an inverse 
filtering whose transfer function would have a form 
approaching the expression 1/R^(Z), where R^(Z) denotes the 
polynomial in Z"^ of order W-p having roots that are the 
W-p roots furthest from the unit circle in the complex 
plane. In particular, it could be of the *^zero- forcing" 
type. This process generates only a small enhancement of 
the noise, since the roots of the associated transfer 
function in Z are relatively far from the unit circle. 

For the p roots which are closest to the unit 
circle, we adopt measures to avoid or to limit the 
incidence of the noise enhancement problem. A MMSE or 
similar method can be chosen as the "second equalization 
method". However, this second method will advantageous 
take into account the discrete nature of the unknown 
symbols. In particular, it may rely on a trellis 
algorithm, such as the Viterbi algorithm, whose 
implementation is common in channel equalizers where the 
system size is not too big. 

The implementation of the second equalization method 
is generally more complex than the first. In each 
particular case, the choice of the number p allows the 
best compromise between the reliability of the 
estimations, which favors higher values of p, and the 
complexity of the equalizer, which favors lower values of 
P- 

Another aspect of the present invention relates to a 
radio communications receiver comprising : 



- conversion means to produce digital samples from a 
radio signal received over a transmission channel 
represented by a finite impulse response of W+1 
coefficients, W being an integer greater than 1; 

- means for measuring the channel impulse response; 

- means for calculating the W roots in the complex 
plane of the Z- transform of the impulse response; 

- means for distributing the W roots into a first set 
of W-p roots and a second set of p roots, p being an 
integer greater than 0 and smaller than W, the roots 
of the second set being closer to the unit circle 
than those of the first set according to a 
determined distance criterion in the complex plane; 

- a first equalization stage for producing an 
intermediate signal by applying a first equalization 
method to the received signal based on a finite 
impulse response whose Z transform, formed by a 
polynomial of degree W-p in Z"-^, has roots which are 
the W-p roots of the first set; and 

- a second equalization stage for producing 
estimations of the discrete symbols of a signal 
carried on the channel by applying a second 
ec[ualization method to the intermediate signal based 
on a finite impulse response whose Z transform 
(R-^(Z)), formed by a polynomial of degree p in Z'^, 
has roots which are the p roots of the second set. 
Other features and advantages of the present 

invention will become apparent in the following 
description, with - reference to the attached figures, of 
non-limiting embodiments, in which : 

- Figure 1 is a block diagram of an example of a radio 
communications receiver according to the invention; 



- Figure 2 is a flow chart showing an embodiment of 
the method according to the invention; and 

- Figure 3 is a diagram illustrating the performances 
of the method- 

The receiver represented in Fig. 1 comprises a radio 
stage 1 which receives the signal picked up by antenna 2 
and converts it to the baseband. The baseband signal is 
digitized by the analog- to-digital converter (ADC) 3 and 
fed to the reception filter 4. Filter 4 provides matched 
filtering with respect to the signal shaping by the 
transmitter. It delivers a digital signal at a rate of one 
complex sample per transmitted symbol. 

The digital signal is fed to a demodulator 
comprising, on the one hand, a synchronization and channel 
estimation module 6 and, on the other hand, an equalizer 
7, 

The synchronization and channel estimation are, for 
example, performed conventionally by means of a 
synchronization sequence included by the transmitter in 
each signal frame. Detection of this sequence, known to 
the receiver, allows, on the one hand, the receiver to be 
synchronized relative to the temporal structure of the 
transmitted frames, and on the other to estimate the 
impulse response r = (rQ, r^^ , . • • , r^) of the channel carrying 

the transmitted frames. The impulse response calculated by 
module 6 is fed to the equalizer 7. 

Equalizer 7 functions for example according to the 
flow chart shown in Fig. 2 to process each frame of the 
received signal, which are presented in the form of a 



vector Y= 



, with L=n+W, adopting the previous notation. 



The channel estimation module 6 having supplied the 



W+1 complex coefficients of the estimated channel 

impulse response, the first step 10 consists of a search 
of the W roots of the Z- transform of the impulse response, 
given by equation (5) . Various conventional methods for 
the determination of complex roots of a polynomial can be 
employed at step 10, To this end, reference can be made to 
the work of E. DURAND : "Numerical Solutions of Algebraic 
Equations, Volume I, Equations of the Type F(x)=0", 
Editions Masson, I960, 

The W roots thus found a^^, a2 / • . w oc^ are then ordered 
to allow distribution into two sets, one containing the 
W-p roots furthest from the unit circle, and the other the 
p roots closest to the unit circle. 

To this end, a distance 5^ is calculated at step 11 
for each of the roots (l<q<W) , This distance can be 

conveniently obtained as follows : 
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>1 



(7) 



In step 12, the roots of the transfer function 

R(Z) are sorted in order of decreasing distance: 
> 82 ^ ... > 5^. Then, the first W-p roots cx-^/ • • • # oc^.pr 
which are the furthest from the unit circle, are separated 
from the remaining roots ct^.p+i/ . . . ,a^. 

In step 13, the equalizer 7 develops a polynomial in 
Z"^ defined by : 

w-p w-p 

R^Z) = n(l-aq-2-^)= Zsa-Z""^ (8) 



q=l q=0 

This allows the coefficients s. 



of the transfer 



function R^(Z) associated with the impulse response 



'W-p 



) of a virtual channel to be determined. 



which corresponds to the estimated transmission channel 
with the contributions closest to the sincfularity zones 
removed . 

We then proceed to a first equalization 14 which is 
effectively an inverse filtering approaching the transfer 
function 1/R^(Z). Several implementations can be used to 
perform this inverse filtering. In particular, a zero- 
forcing equalization can be performed, as previously 
discussed. Regarding these methods, reference can be made 
to the work of J.G. Proakis : ^^Digital Communications", 
McGraw-Hill, 2^^ Edition, 1989. 

The inverse filtering 14 produces an intermediate 
signal in the form of a vector Y' of L'=n+p samples 
y' • • • 'Y'l* case of the zero-forcing method, the 

vector Y' is obtained from the matrix equation : 



Y' = (A'^ A' ) A'" Y 



(9) 



In expression (9), A' denotes a matrix having a 
Toeplitz structure with n+W rows and n+p columns, 
generated from the coefficients s^ of the polynomial R^(Z): 



A' = 



0 
0 
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the roots 



(10) 



results in the 



eigenvalues of matrix A'^A' being relatively far removed 
from 0 . 

Alternatively, the inverse filtering can be achieved 
by cascading W-p filtering cells, each corresponding to 



-li- 
the inverse of a transfer function (Z) = 1 - a Z'^, for 



l<q<W-p. If I Ug I =1, the inverse filter of (Z) is 

unfeasible. If | | < 1, we can develop 1/R^ (Z) in the 
form : 

— = 1 + a^, z"^ + a^. Z"^ + . . . + a^- Z""" + . . . (11) 

R|(Z) q q q 

Development (11) is causal and stable since the 
convergence domain contains the unit circle. The inverse 
filter cell can thus be generated in either a transverse 
or a recursive form. 



If I ttq I > 1/ then 1/R^ (Z) can be developed in the 
form : ^ 



1 -1 / -1 

= = -Oifj . Z . 1 + ttcr 



s - ^„ • " ^ • Z + a^^. Z^ + . . . -h aj^. + . . .) (12) 
Rq(Z) ^ V ^ ^ ^ 

Development (12) is anti -causal and stable. To perform the 
inverse filtering cell, the development (12) is truncated 
and an transverse filter implementation is adopted. The 
anti -causality results in a delay corresponding to the 
length of the retained response. 

We note that developments (11) and (12) justify the 
unit circle distance criterion 5^ used according to 
relation (7) . 

In step 15/ the equalizer 7 develops a polynomial of 
degree p in Z"-^, whose roots correspond to the p roots of 
R(Z) closest to the unit circle, such that 
R(Z) = rS(Z) .R^(Z) : 

R^(Z) = n - Otq- Z = 2 tq- Z""^ (13) 

q = W-pH-l q=0 

The complex coefficients t^ define the impulse 
response of another virtual transmission channel, whose 
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equalization by a zero- forcing method or the like would 
pose problems of noise enhancement. 

The intermediate signal Y' is then subjected to an 
equalization according to another method, based on the 
impulse response t = {t^, t^, . . . , t^) . This second 
equalization 16 can be conveniently performed by means of 
a Viterbi trellis (see above-cited article from G.D. 
Forney Jr., or the above-cited work of J.G. Proakis) . 

The second equalization stage 16 generates 
estimations dj^ of the symbols of the frame (l<k<n) . These 
estimations dj^ constituting the output of equalizer 7 can 
be input to a deinterleaving module 8 then to a channel 
decoder 9 which detects and/or corrects possible 
transmission errors . 

Figure 3 illustrates the performance of the method 
in the case of the transmission of a signal frame 
according to the GSM European mobile telephone system, 
replacing the binary modulation of the GMSK type by an 
eight state phase modulation (8-PSK modulation) . The 
impulse response was truncated at five bit time periods 
(W=4) . Figure 3 shows the dependence of the bit error rate 
(BER) , expressed in %, on the ratio Sj^/Nq between the 

energy-per-bit and the spectral density of the noise, 
expressed in decibels (dB) • The BER is that observed in 
the symbol estimations after deinterleaving and channel 
decoding carried out by a method conforming to those 
employed in the GSM standard. Curve I shows the results 
obtained by the pure zero -forcing method, i.e. in the 
limit case where p=0. Curve II shows the theoretical 
result that would be obtained by equalizing the channel 
purely with the Viterbi algorithm (limit case where p=W) • 
In practice, the corresponding trellis should have 8^=4096 
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states, meaning that the Viterbi equalizer could not - be 
realized using current technology. The difference between 
curves I and II illustrates the superiority of the Viterbi 
algorithm which produces maximum likelihood estimations. 

Curves III and IV show the results obtained by the 
method according to the invention, in the cases where p=l 
and p=2, respectively. The quite significant improvement 
already achieved for the value p=l with respect to pure 
zero forcing can be seen. 

By way of an example, equalization of a GSM signal 
frame by a pure Viterbi algorithm, under the conditions of 
Figure 3, would require of the order of 8.45 million 
floating-point operations, i.e. around 1.83 Gflops, while 
the implementation of the present invention under the same 
conditions requires of the order of 19,000 floating-point 
operations («4 . 2 Mflops) in the case of p=l, including 
searching the roots of R(Z) and 1/R^(Z) inverse filtering 
by the zero- forcing method. The niunber increases to the 
order of 129, 000 operations («28 Mflops) in the case of 
p=2, which remains compatible with the capacity of 
currently available digital signal processors (DSP) . 



